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We present a density functional theory study of low density bromination of graphene and 
graphite, finding significantly different behaviour in these two materials. On graphene we find 
a new Br2 form where the molecule sits perpendicular to the graphene sheet with an extremely 
strong molecular dipole. The resultant Br"'"-Br~ has an empty p^-orbital located in the graphene 
electronic vr-cloud. Bromination opens a small (86meV) band gap and strongly dopes the graphene. 
In contrast, in graphite we find Br2 is most stable parallel to the carbon layers with a slightly 
weaker associated charge transfer and no molecular dipole. We identify a minimum stable Br2 
concentration in graphite, finding low density bromination to be endothermic. Graphene may be a 
useful substrate for stabilising normally unstable transient molecular states. 

PACS numbers: 31.15.A- 73.22.Pr 81.05.uf 37.30.+i 



The intercalation of carbon nanomaterials with elec- 
tron donors and acceptors is an active research area 
in which much effort is channelled towards the under- 
standing and controlling of the electronic properties of 
graphene and graphite. Numerous potential applications 
such as sensors, electronic display panels, hydrogen stor- 
age and supercapacitord^ have been suggested for such 
intercalated materials. The layered structure of graphites 
plays an important role in charge transfer reactions. Ac- 
ceptor species can intercalate between graphitic layers, 
expanding the graphite with the resultant hybrids known 
as graphite intercalated compound^. 

Bromine acts as an acceptor when intercalated in ma- 
terials such as graphite or nanotubes, and has been pro- 
posed experimentally as a way to open a band gap in 3- 
or 4- layer graphene^l. The in- plane electrical conduc- 
tivity of graphite increases from 2.4 x 10'' ^l~^cm~^ at 
room temperature to 2.2 x 10^ f2~^cm^^ after intercala- 
tion with bromin^. 

Bromine forms many ordered phases in graphites 
and undergoes an order-disorder phase transition as 
the amount of bromine or temperature changeJ^^^^. 
Bromine intercalated graphite forms stage-n compounds 
where n is the number of graphitic layers between planes 
of Br2 {n > ^ . Extended X-ray Absorption Fine Struc- 
ture (EXAFS) spectroscopy at room temperature and 
lOOMi^ showed that intercalated bromine molecules lie 
parallel to the basal plane, with an expansion of the Br- 
Br distance by 0.03 A to accommodate the lattice mis- 
match between the free molecule and the 2.46 A spac- 
ing between graphite hexagons. X-ray diffraction and 
electron microscopy studiegSl suggest that intercalated 
Br2 at lower concentration is composed of chains of Br2 
molecules in which the intermolecular distances is iden- 
tical to that of solid bromine. 



Graphite Raman spectra associated with Br2 intercala- 
tion show a strong peak at 242 cm~^Ii^lti^ assigned to the 
intercalated Br2 stretch mode. The frequency is down- 
shifted from 320 cm~^ for gaseous Br2 and 295 cm~^ 
for solid molecular bromine^. There have been limited 
density funct ional studies of brominated graphitd^ and 
graphen d^ l ^^ l 

In this paper we examine low density bromination of 
graphene and graphite using density functional (DFT) 
calculations within the local density approximation^-^. 
The method has been successfully used to study inter- 
calated boron in graphitd^. A localised Gaussian basis 
set is used with a large number of fitting functions per 
atom (22 for each C atom and 50 for each Br), with 
angular momenta up to 1=2 for C and 1=3 for Br. A 
finite temperature electron level filling of kT=0.04eV is 
used to improve convergence. Core electrons were elim- 
inated using norm-conserving relativistic pseudopoten- 
tials of Hartwigsen, Goedecker and Huttep^Sl. A cut-off 
energy of 150 Hartrees was used to obtain convergence of 
the charge density. 

Isolated Br2 was calculated in a 13.23 A cubic su- 
percell. Hexagonal 4x4 graphene supercells containing 
C32Br2 were used with a large vacuum spacing of 31 A 
between layers to ensure no inter- layer interaction, and 
a 4 X 4 X 1 Monkhorst-Pack k-point gricP^^. Graphite cal- 
culations used 3 X 3 X n layer supercells (Ci8)„(Br2)m, 
n=l-4, 171=1-2) for different layer stackings, with 4x4x1 
or 4 X 4 X 2 k-point grids depending on cell size. All 
structures were fully geometrically optimised with no 
constraints of symmetry, allowing both atomic positions 
and cell dimensions to vary freely. Atomic charge states 
were obtained by summing MuUiken population analysis 
over all the filled electronic states. Vibrational frequen- 
cies were calculated by determining the energy and forces 
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FIG. 1. Optimised geometries of Br2 molecule perpendicular 
(a-c) and parallel (d-h) to a graphene sheet. Relative stabili- 
ties are quoted in eV (binding energy of Br2 in structure (a) 
is 0.40eV). The text refers to (a) Br^ and (d) Br!| 



TABLE I. Calculated parameters for the most stable Br2 and 
Brj on graphene compared with literature DFT calculations. 





Br^ 


Br| 


Literature 


Orientation 


Perp. 


Parallel 


Parallel 


Binding energy / Br2 (eV) 


-0.40 


-0.29 


-0.2dil 








-0.2<P 


Br-Br Stretch Frequency (cm~^) 


288 


270 


3110 


Br-Br bond length (A) 


2.33 


2.31 


2.24^ 


Br-C distance (A) 


2.74 


3.45 


3.37^ 








3.7#3 


Charge state of Br atoms (e) 


+0.48 


-0.04 






-0.61 


-0.04 




Total Charge transfer / Br2 (e) 


-0.13 


-0.08 





for ± 0.2 au displacements of the Bromine atoms. The 
second derivatives on the displaced atoms can then be 
found from the two-sided difference formula for the sec- 
ond derivative. All results are spin averaged, test calcu- 
lations with spin polarisation all gave zero spin solutions 
as the most thermodynamically stable. 

Our calculated bond length, 2.29A, and stretching 
frequency, 326 cm~^, for isolated Br2 show excellent 
agreement with experiment (2.27A and 323 cm.~'^^, 
2.283A and 320 cm"^ and literature DFT/LDA values 
(2.263AI^, 2.244A and 324 cm"^. 

Figure [T] shows the structures of Br2 over graphene in 
both perpendicular and parallel orientations after geom- 
etry optimisation, with calculated properties for Brn_and 

Br| (structures (a) and (d) respectively in Figure llj) in 
Table H 

The most thermodynamically stable arrangement is 
Br2 oriented perpendicular to the graphene sheet above 
a carbon atom (Figure Im, referred to hereafter as Br^ 



). Its binding energy of 0.40eV shows it will be strongly 
physisorbed at room temperature. However the small 
variations in binding energy between structures suggests 
that low density Br2 binding to graphene will be largely 
orientation independent; indeed and at these densities at 
room temperature Br2 should be in constant tumbling 
motion (this also holds for the results oP^. Increasing 
Br2 concentration did not significantly change the rela- 
tive energies of perpendicular and parallel orientations. 
However it is possible to obtain twice the maximum sur- 
face density for Br^ as for Br2 in the parallel orientation 

(Figure [iji, referred to hereafter as Br2 ). Thus for the 
limit of high surface concentrations we expect Br^ to 
dominate. 

Br^ represents a very unusual configuration for 
bromine. It shows strong charge transfer (0.129e) from 
the graphene, with a very strong induced molecular 
dipole (Br+'' ''8°-Br-°-6°9). The singly occupied p,- 
orbital of the lower Br atom depopulates into the Pz of 
the upper Br atom forming a nascent Bromonium and 
Bromide ion pair. In this way the emptied lower p^- 
orbital can sit within the graphene 7r-cloud (Br only 2.74 
A above the graphene). 

This behaviour is reminiscent of the well-known reac- 
tion between Br2 and unsaturated bonds in organic chem- 
istry. However in these cases this dipolar form of Br2 is an 
unstable transient state that immediately saturates the 
C=C bond, forming two Br-C bonds. In graphene this 
final step would be endothermic due to steric hindrance 
between the Br atoms as a result of the mechanical con- 
finement of the lattice. Indeed our attempts to stabilise 
C-Br pairs on graphene in both neighbouring (1,2) and 
cross-hexagon (1,4) configurations both resulted in Br 
spontaneously reconstructing into a Br2 molecule (plac- 
ing Br-C in a (1,4) configuration with Br atoms on oppo- 
site sides of the graphene is 1.68eV less stable than Br^ 
)• 

The unusual Br^ configuration is reflected in the band 
structure (Figure [2^). The strong coupling between the 
graphene LUMO and the Br2 antibonded state at around 
+0.6 eV reflects the interaction between the empty Br p^- 
orbital and the graphene 7r-cloud, with the resultant low 
density Br2 layer opening a small 86meV band gap. 

In contrast Br| has no induced dipole, with weaker 
charge transfer (0.084e) from the graphene. The molecule 
sits 3.254 A above the graphene i.e. above the 7r-cloud, 
with the additional charge occupying the Br2 ppa* anti- 
bonding state. The band structure (Figure [2]d) shows the 
Br2 states lie lower than those of Br^ by ~0.7eV. The 
Br2 anti-bonded state pins the Fermi level '-^0.2eV lower 
than in the pristine case, leaving graphene states around 
the K-point depopulated, indicating charge transfer from 
graphene to Br2 . The Bromine states are flat and largely 
decoupled from the graphene bands since there is only 
weak interaction between Br2 orbitals and the graphene 
TT-cloud in this orientation. 

Our results are consistent with experimental Ra- 
man observations of brominated graphen^. While for 
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FIG. 2. Band structure of (a) Br^ and (b) Brj on graphene 
(eV) , faded dotted lines indicate the same supercell of pristine 
graphene for comparison. Arrows indicate bromine related 
states. 



graphene with three or more layers a resonant Raman sig- 
nal was observed for Br2 at around 240cm~^, for mono- 
and bi-layer graphene no bromine signal was observed. 
This could be an orientation effect, since if the Br2 sits 
perpendicular to the surface as we propose and orthog- 
onal Raman is used, then the molecules will be aligned 
with the beam and there will be no interaction and hence 
no signal. In addition the bromine HOMO/LUMO states 

are further apart for Br^ (^2.2eV) than for Br2 (^1.96) 
or the various graphite structures we have examined (1.7- 
2.0eV). Given the excitation laser used (633nm, 1.96eV), 
Br^ may not be in resonance as the authors suggestecP. 
Br| would have an associated Raman signal, and since 
none is observed this allows us to exclude this configura- 
tion. We note that the strong dipole of Br^ will make 
the molecule infra-red active. 

We now turn our attention to graphite. Our cal- 
culated energy to separate AB graphite layers is 36.7 
meV/atom, with AA stacked graphite 12.0meV/atom 
less stable than AB stacked graphite. These figures are in 
good agreement with experiment (35meV/atonP2l) and 
previous calc ulation s (9.68-9.70meV/atom AA/AB en- 
ergy difFerenc^2ll2ll) _ Our interlayer spacings of 3.39A 
and 3.50A for AB and AA-stacked graphite respec- 
tively are also in reasonable agreement with previous 



TABLE II. Calculated and experimental results for Br2 inter- 
calated graphite. The C-C layer distance refers to the layers 
separa ted by bromine . Exper imental values from EXAF^, 
Ramarji3Iir and XREPMIMI. At these low concentrations 



Br2 intercalation is endothermic. 





Stage- 1 


Stage-2 


Experiment 


Cell used 


CisBrz 


C72Br4 




Binding Energy /Br2 (eV) 


-1-0.08 


+0.08 




C-C layer distance (A) 


6.45 


6.47 


7.d^7.0d^ 


c-axis (A) 


9.80 


9.82 


10.#> 


Br-Br bond length (A) 


2.30 


2.30 


n 04IOI27II2I 


Br-C distance (A) 


3.51 


3.46 


2.cP 


Br-Br frequency (cm~^) 


287 


274 


242-25^^^1^^ 


Charge transfer / Br2 (e) 


-0.10 


-0.12 


-0.1tP-0.34Sl 



calculation^2ll. 

We placed Br2 in a variety of different orientations and 
locations including above a- and /?- carbon atom sites 
and hexagon centers, in graphite of various layer stack- 
ings. Unlike graphene, Br2 in graphite is more stable par- 
allel to the graphitic layers with Br atoms above hexagon 
centers (see Table |ll]), in agreement with experimen|S-i21. 
Perpendicular oriented Br2 structures are much less sta- 
ble (by typically 0.52eV). The results are summarised 
and compared with experiment in Table[ll] We note that 
the charge transfer value from Raman was reported with 
a large uncertaintjff^. 

Besides the lowest energy structures quoted here we 
obtained many metastable structures. While their ener- 
gies were all within 0.01-0.05eV of those structures dis- 
cussed here, they show significant variation in the Br2 
stretch frequency (250-278cm~^), and slight variation in 
position and Br-Br bond length (2.31-2.33A). This sug- 
gests that at room temperature Br2 in graphite will be 
mobile, and is consistent with the observation of a broad 
and somewhat complicated Raman pealsP^. 

We explored all possible layer stacking combinations 
for stage-1 and stage-2 intercalated graphites. Bromine 
molecules are most stable with AA-stacked graphite each 
side, while unbrominated graphite layers preferentially 
are AB stacked. Thus (indicating layers of bromine 
molecules by X) we find the most stable stage-1 phase 
to be [AX]„, and stage-2 to be [AXABXB]„, suggesting 
[AXAB]„ and [AXABABXBAB]„ for stagc-3 and stage-4 
respectively. 

Our calculated bromine intercalation energy is weakly 
endothermic, since at these low densities the energetic 
cost associated with separating graphite planes (a cost 
per unit area) is not sufficiently offset by the binding 
energy of Br2 to the layers. Subtracting the energy to 
separate graphite layers from the Br2 intercalation en- 
ergy gives an energy for intercalation of Br2 into 'pre- 
separated stage-1 graphite' of 0.581eV/Br2. This implies 
that the minimum Br2 concentration for exothermic in- 
tercalation in stage-n graphite will be Ci6nBr2. Indeed 
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a fixed-Br calculation for a C8Br2 stage-1 high coverage 
structure gives Br2 intercalation as exothermic. Thus Br2 
will aggregate in the same inter-layer space and should 
be considered as a layer rather than individual molecules. 

This minimum required concentration also suggests 
intercalation will be a slow diffusion process with an 
abrupt diffusion front. This is consistent with long exper- 
imental intercalation times^*, despite the high bromine 
inter-layer mobility. It also explains why, on out-gassing 
bromine, the material switches from a stage-2 to com- 
pound to stage-n (n=3, 4, ...) rather than remaining 
stage-2 with lower bromine density per layeJ^EHSll. 

We note that frequency calculations incorporating the 
energetic double derivatives of surrounding carbon atoms 
gave identical values to within 1 cm~^, showing that the 
Br2 stretch mode is decoupled from the surrounding car- 
bon lattice consistent with the literatur^. 

Band structure calculations of Br2 layers between 
graphite sheets (not shown here) give a bromine-related 
state which pins the Fermi level just below that of per- 
fect graphite (~0.1eV) indicating charge transfer from 
graphite to Br2. This state shows some dispersion indi- 
cating weak coupling with the underlying graphite. In 
other respects the graphite band structures are barely 
perturbed. 

Our results on graphite and graphene can explain the 
anomalously large G peak shift for single graphene sheets 
in comparison with few-layer graphene^J. Since the higher 
maximum Br2 surface density for Br^ than for inter-layer 
Br2 means charge transfer per unit area will be higher. 
Additionally Br2 can attach to both sides of graphene. 
We find a binding energy of -0.38eV /Br2 for two Br^ 
either side of the same C atom, with associated charge 
transfer of -0.12e/Br2. Thus the net total charge trans- 
fer per unit area will indeed be significantly higher for 
monolayer graphene than multi-layer systems. 



Test calculations for a (5,5) single walled nanotube, 
either isolated or in bundles, gave similar structural be- 
haviour, i.e. Bi2 on the surface of the isolated tube 
adopts a perpendicular orientation, while intercalated 
Br2 sits parallel to the tube walls. This will be explored 
further in a later publication. 

In summary, we have examined low density Br2 adsorp- 
tion in graphene and graphite. On graphene Br2 adopts 
an unusual perpendicular orientation, opening a small 
band gap (^86meV) in the graphene with strong charge 
transfer. The molecule forms a Br^-fBr^ pair, rendering 
it infra-red active. This is a new form of Br2 previously 
only considered as an unstable intermediate to bromine 
induced carbon bond saturation. Such graphene-induced 
stabilisation behaviour may be mirrored in other molecu- 
lar species, enabling study of otherwise unstable reactive 
molecular forms. 

In graphite Br2 adopts a parallel orientation to the 
sheets with an associated charge transfer. Our calcu- 
lations are in good agreement with experimental data 
where available. At high bromine concentrations and 
low temperatures t here i s some evidence of bromine chain 
structure format iorP^ESlin graphite, and we are currently 
investigating this further. We note that high density 
bromination of graphite leads to stage-2 compounds, and 
in conjunction with an appropriate secondary surfactant 
this may be a promising way to produce bilayer graphene. 
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